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Aims. Physical and evolutionary properties of the sub-mJy radio population are not entirely known. The radio/optical analysis of 
the ATESP 5 GHz sample has revealed a significant class of compact flat/inverted radio- spectrum sources associated to early-type 
galaxies up to redshift 2. Such sources are most plausibly triggered by an AGN, but their observational properties are not entirely 
consistent with those of standard radio galaxy populations. In the present work we aim at a better understanding of the radio spectra 
of such sources and ultimately of the nature of AGNs at sub-mJy flux levels. In particular we are interested in assessing whether the 
(~ l , AGN component of the sub-mJy population is more related to efficiently accreting systems - like radio-intermediate/quiet quasars - 

^ I' or to systems with low accretion rates - like e.g. FRI radio galaxies - or to low radiative efficiency accretion flows - like e.g. ADAF. 

I ' Methods. We used the ATCA to get multi-frequency (4.8, 8.6 and 19 GHz) quasi-simultaneous observations for a representative 

O ' sub-sample of ATESP radio sources associated with early-type galaxies (26 objects with 5 > 0.6 mjy). This can give us insight into 

the accretion/radiative mechanism that is at work, since different regimes display different spectral signatures in the radio domain. 
^ . Results. From the analysis of the radio spectra, we find that our sources are most probably jet-dominated systems. ADAF models are 

^ I ' ruled out by the high frequency data, while ADAF+jet scenarios are still consistent with flat/moderately inverted-spectrum sources, 

but are not required to explain the data. We compared our sample with high ( ; 20 GHz) frequency selected surveys, finding spectral 
properties very similar to the ones of much brighter (5 > 500 mJy) radio galaxies extracted from the Massardi et al. (2008) sample. 
^ ■ Linear sizes of ATESP 5 GHz sources associated with early type galaxies are also often consistent with the ones of brighter B2 and 3C 

' radio galaxies, with possibly a very compact component that could be associated at least in part to (obscured) radio-quiet quasar-like 

' objects and/or low power BL Lacs. 

^ : Key words. Surveys - Rad>o contmuum: general - Methods: data analysis - Catalogues - Galax.es: general - Galaxies: evolution 



1 . Introduction intermediate/quiet quasars - or to systems with very low accre- 

— . . , tion rates - like e.g. FRI radio galaxies (FanarofF & Riley I97J). 

O After many ye^s of studies, AGNs are now recognised to ^j^^ j^^^^^ ^^^^^^.^ ^^^^.^ ^^^^^ supported by the presence of 

contribute significantly at radio fluxes below 1 milliJy (mJy). j^^jj .^^^^^^ j^^^^^ ^ 

Multi-wavelength studies of deep radio fields show that star- ^^^.^ ^^^^^^^ ^^^^ ^^^^^^.^^ ^^^^^ ^^^^ 

. ^ . N ' ^Jy^ l^^^lf (Richards et .^^pp^ncd by the large number of so-called radio-intermediate 

H . al. 1 1999) while radio sources associated with early-type galax- ^^^^^ ^^^^^.^^^ j^^^j^ ^^^^ Lacy etal. IMTl) and by 

les and plausibly triggered by AGNs, are the most significant ^j^^ modelling work of Jarvis & Rawlings (&, who predict a 

source component (- 60 - 70% of the total) at higher flux significant contribution of radio-quiet quasars at sub-mJy levels, 

densities (> 100 - 200 A-Jy) with a further 10% contribution ^^^^^j^^^ .^^^^ represented by the possible role played at low 

from broad-/narrow-line AGNs (see e.g. Gruppiom et al. d999t ^^^-^ ^^^^^ ^ ^^^-^^-^^ efficiency accretion mechanisms, 

Georgakakis al. 1999; Magliocchetti et a . 2000, Prandoni et ^,,ociated to optically thin discs, such as the so-called advection 

" • i^^^i o f " ? .^^o^' ;i f ™^ dominated accretion flows (ADAF) and modifications (ADIOS, 

al. \2m; Smolcic 2008). This somehow unexpected pres- ^.j^^p^ ^ YiCM Quataert & NarayanflM 

ence of large numbers of AGN-type sources at the sub-mJy level Abramowitcz et al 206^ 

has given a new and interesting scientific perspective to the study ^ ■„•,,■ , , , , 

of deep radio fields, since a better understanding of the physi- , , ^^^^ especiaUy suited radio sample to study the phenomenon 

cal and evolutionary properties of such low/intermediate power °f low-luminosity nuclear activity, possibly related to low radia- 

AGNs may have important implications for the determination of tion/accretion processes and/or radio-intermediate/quiet QSOs, 

the black-hole-accretion history of the Universe as derived from ^^ESP 5 GHz survey carried out with the 

radio-selected samples. Of particular interest is the possibility of ^^^'^ Telescope Compact Array (ATCA) by Prandoni et al. 

assessing whether the AGN component of the sub-mJy popula- <™, which coyei^ a 2 x 0.5 sq. degr region^e wi^r 

tion is more related to efficiently accreting systems - like radio- ^^^'lu c^Jf.o'"^'^ ^'^'"'"''""^ et ai ^^O^- 
2001a). The ATESP 5 GHz survey consists of 131 radio sources 

Send offprint requests to: I. Prandoni, e-mail: with flux densities larger than S ~ 0.4 mJy detected either at 

i.prandoni@ira.inaf.it 1.4 GHz or 5 GHz or both. Interestingly, the analysis of the 
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1.4 and 5 GHz ATESP surveys has revealed a significant flat- 
tening of the source radio spectra going from mJy to sub-mJy 
flux levels (Prandoni et al. 2006). Such flattening is mostly asso- 
ciated to early-type galaxies as demonstrated by deep {R < 25) 
UB VRIJK multi-colour imaging. At the flux limit of the sample 
star-forming galaxies are starting to appear, but are not yet the 
dominant population: ~ 14% of the sources are identified with 
broad/narrow-line AGNs and another ~ 64% is identified with 
early-type gal axies, most probably triggered by an AGN (see 
Mignano et al. l2008h . The absence of emission lines in the latter 
together with low radio luminosities (typically 10^^"^^ W/Hz) 
suggests that these are FRI radio galaxies. However, > 60% 
of them show flat (a > -0.5, where S ~ v^) and/or inverted 
(a > 0) spectra and rather compact linear sizes (d < 10-30 
kpc, see Mignano et al. 2008). Both compactness and spectral 
shape suggest a core emission with strong synchrotron or free- 
free self-absorption. Such sources are known to exist among FRI 
radio galaxies, but they are relatively rare. It is therefore impor- 
tant to confirm the existence of such a large class of flat/inverted 
spectrum objects through simultaneous radio observations at dif- 
ferent radio wavelengths. 

If real, these sources may represent a composite class of ob- 
jects very similar to the so-called low power iP4omHz < 10^^^ 
W/Hz) compact (< 10 kpc) - LPC - radio sources studied by 
Giroletti et al. (l2005b . LPC host galaxies do not show signa- 
tures of strong nuclear activity in the optical (and X-ray) bands, 
and preliminary results indicate that multiple causes can pro- 
duce LPC sources: geometrical-relativistic effects (low power 
BL-Lacertae objects), youth (GPS-like sources), instabilities in 
the jets, frustration by a denser than average ISM and a prema- 
ture end of nuclear activity (sources characterised by low ac- 
cretion/radiative efficiency, i.e. ADAF/ADIOS systems). Some 
difficulties remain however. First, if the objects are young GPS 
sources, much higher (> 10^^ W/Hz) radio luminosities are ex- 
pected; second, if they are old ADAF/ADIOS sources, the ex- 
pected luminosities are much lower (< 10^' W/Hz, Doi et al. 
I2005h . In the latter case there may yet be different solutions, as 
for example if the ADAF source coexists with a radio jet, as was 
suggested by Doi et al. (2005) in the case of a few low luminos- 
ity A GNs (see also the ADAF-jet model of Falcke & Biermann 
1199^) . 

Another possibility is that we are dealing with radio- 
quiet/intermediate QSOs, where the activity in the optical band 
is obscured by dust in the galaxy. This scenario may be sup- 
ported by the results reported by Klockner et al. (2009), 
who have recently observed with the EVN 11 z > 2 radio- 
intermediate obscured quasars (several with flat-spectrum), de- 
tecting seven of them. The detected radio emission accounts for 
30-100% of the entire source flux density, and the physical ex- 
tent of this emission is ; 150 pc. The missing flux implies the 
likely existence of radio jets of physical sizes between ! 150 pc 
and ;40 kpc. 

The aim of the present paper is a better understanding of the 
physical properties of the AGNs associated to early-type galax- 
ies in the ATESP 5 GHz sample. The study of the radio spec- 
tra, based on multi-wavelength (simultaneous) observations, is 
essential in this respect, since different source types (e.g. GPS 
vs BL Lac) and/or accretion scenarios (e.g. jet-dominated vs 
ADAF/ADIOS -dominated) lead to different shapes of the radio 
spectrum. Particularly useful is the availability of high radio fre- 
quency information, since, for example the location of the fre- 
quency at which the radio spectrum peaks may be an indication 
of the presence or absence of outflows. A pure ADAF model 
without outflow is expected to peak in the sub-mm (the "sub- 



mm bump") and has an inv erted (a > 0) spectrum in the cm 
range (see e.g. Nagar et al. 2001*). As outflows become impor- 
tant, the peak shifts to longer wavelengths (above 1 cm), and de- 
pending on which frequencies were observed, the spectra might 
be classified as flat [a > -0.5) or steep (a < -0.5). 

In this paper we report on quasi-simultaneous observations 
at 4.8, 8.6 and 19 GHz of a subset of ATESP 5 GHz sources asso- 
ciated to early-type galaxies. Our results will be compared to the 
ones obtained by new high frequency (> 10 GHz) radio surveys. 
Such surveys have starte d only recently, nota bly wi th the ATCA 
at 20 GHz (Ricci et al. 120041: Sadler et al. 120061: Massardi et 
al. [2008), and detailed spectral studies based on high radio fre- 
quency observations are therefore still scarce. Other dee per re - 
cent high frequency surveys are the one by Tucci et al. (l2008l) . 
who observed sources with 5 > 20 mJy at 33 GHz and compared 
the composition of the source counts with the high frequency 
model counts of De Zotti et al. (2005), and the 9C survey at 15 
GHz (a 124 source sample complete to 25 mJy an d a 70 source 
sample complete to 60 mJy, Bolton et al. |2004|) . In the latter 
case the observed source counts were discussed in Waldram et 
al. (2007). The data presented in this paper represent the first 
high frequency (> 10 GHz) systematic follow-up of a sample of 
sub-mJy radio sources. 

The paper is organized as follows. We discuss the selection 
of the sample, the follow-up observations at 4.8, 8.6 and 19 GHz 
and the data reduction in Sects. |2] and [3l source variability at 
4.8 GHz (where data at two epochs are available) is discussed in 
Sect.m a complete discussion of the radio spectra, in comparison 
with high-frequency selected samples, is given in Sect. |5] A fur- 
ther comparison with (FRI and FRII) radio galaxies and radio- 
quiet QSOs is presented in Sect. |6l where we discuss source lin- 
ear sizes and radio powers. A brief summary of the main results 
of this work is given in Sect. IT] Throughout this paper we use 
the ACDM model, with Hq = 70, Q„, = 0.3 and Q.a = 0.7. 



2. Sample selection 

In order to better assess the composite nature of the sub-mJy 
population with particular respect to the low-/intermediate- 
luminosity AGN component, we are studying in great detail a 
sample of 131 sources with5 > 0.4 mJy distributed over a 2x0.5 
square degree region in the Southern sky surveyed at both 1 .4 
and 5 GHz in the framewor k of the ATESP 1.4 an d 5 GHz ra- 
dio surveys (Prandoni et al. l2000al: l2000bl: 1200 laL 12006). This 
region overlaps entirely with the so-called 'Deepl' sub-region 
of the ESO Deep Publi c Survey (DPS) survey (Mignano et al. 
2008; Olsen et al. 2006). As described in detail in Mignano et 
al. (2008), deep (R < 25), multi-colour (UB VRIJK) optical/NIR 
data have been collected for a part of this region (1 xO. 5 sq. degr., 
namely fields Deep la, and b) in the framework of the DPS and 
for another part (0.5 x 0.5 sq. degr., namely field Deeplc) in 
the framework of the Garching-Bonn Deep Survey (GaBoDS, 
Hildebrandt et al. ,2006 ), allowing us to derive photometric red- 
shifts and types for virtually all the radio sources opticall y iden - 
tified down to / < 23.5 in this region (Mignano et al. |2008|) . 
Such data are complemented in the remaining 0.5 x 0.5 sq. degr. 
field (Deepld) by shallower optical imaging (/ < 22.5) obtained 
in th e fram ework of the ESO public survey EIS-WIDE (Nonino 
et g/. 1 19991) and shallower spectroscopic infor mation ( b , < 19.4, 
Vettolani et al. IT998I: / < 19, Prandoni et al. 120011^ : / < 21.5 
Prandoni et al. in prep). This means that in field Deepld opti- 
cal identification and spectral classification is available only for 
~ 33% of the radio sources. 
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Table 1. ATESP-Deepl sample: statistics. Table 2. Follow-up source list observed in May 2007. 
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a,b,c 


85 


56 


37 


48 


32 


23 


23 3*"' 
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46 


15 


10 


31 


9 


5 


3(.) . 


All 


131 


71 


47 


79 


41 


28 


26 3 



not been confirmed 
<*> J225048-400147; J225321-402317; J225504-400154 
J2245 16-40 1807; J224547-400324; J224654-400107 



For the purposes of this work we are interested in the sources 
identified with early type galaxies. As reported in Table [T] 
(Cols. (2)-(4)), we have 85 radio sources in fields Deepla,b,c, 56 
of which have been optically identified and spectrally classified 
(/ < 23.5); 37 of these sources (27 having flat/inverted spectra) 
are associated to early-type galaxies. To these sources we can 
add another 15 (out of 46) sources optically identified and spec- 
trally classified in field Deep Id, ten of which are associated to 
early type galaxies (five with flat/inverted spectrum). 

Subsequently we will focus our analysis on a bright sub- 
sample: for sensitivity reasons we limited our multi-frequency 
(4.8, 8.6 and 19 GHz) follow-up to radio sources with 
5(5GHz) > 0.6 mJy. We have a total of 79 such sources (48 
in fields Deepla,b,c plus 31 in field Deep Id), 41 (32 -i- 9) having 
an optical spectral classification and 28 (23 -i- 5) being classi- 
fied as early-type (see Table[Tl Cols. (5)-(7)). For the 6/3 cm and 
12 mm follow-up (see Table [1] Cols. (8)-(9)) we gave priority 
to the 23 sources in fields Deepla,b,c, whose optical identifica- 
tion is virtually complete down to / ~ 23.5. Nevertheless multi- 
frequency (6/3 cm and 12 mm) quasi-simultaneuous observa- 
tions have been undertaken also for a sub-set of three (out of five) 
Deep Id sources and for three radio sources in fields Deepla,b,c, 
which were preliminarly classified as early-type at the time of 
the 2007 observations, but not confirmed as such in the final op- 
tical analysis (see Mignano et aZ. l2008l) . 

We notice that no target selection was done based upon the 
source radio spectral shape. In other words, both flat/inverted 
and steep spectrum sources were observed. A comparison be- 
tween flat, inverted and steep sources is important to either as- 
sess the possibly composite AGN component of the faint radio 
population, or to test for various accretion models, from pure 
ADAF to ADAF+outflows/jets to pure jet-related accretion. 

We also notice that repeating the 5 GHz observations (al- 
ready available from the ATESP 5 GHz survey) allows us to test 
the importance of variability in the present sample (see Sect.|4|. 

3. Radio observations and data reduction 

3.1. Observations 

The 4.8, 8.6 and 19 GHz observations were carried out in two 
separate campaigns in 2007 (15 sources) and 2008 (14 sources) 
with the ATCA. The expected flux densities at 8.6 and 19 GHz 
were computed based upon the source flux densities at 5 GHz 
and the spectral index between 1.4 and 5 GHz (obtained from 
the original ATESP 1.4 and 5 GHz surveys, Prandoni et al. 
[2000b; 2006). From the expected flux densities, sensitivity and 
total time to be spent on source were derived by requesting at 
least a S/N = 6. 



source name 




R. 


A. 




Dec 


tc/x 


tK 




(12000) 


(12000) 


(min) 


(min) 


J225034-401936 


22 


50 


34.61 


-40 


19 


36.3 


5.0 


4.0 


J225048-400147 


22 


50 


48.04 


-40 


01 


47.0 


22.5 


10.0 


J225057-401522 


22 


50 


57.79 


-40 


15 


22.5 


7.5 


6.0 


J225223-401841 


22 


52 


23.82 


-40 


18 


41.9 


47.5 


42.0 


J225249-401256 


22 


52 


49.92 


-40 


12 


56.0 


35.0 


40.0 


J22532 1-4023 17 


22 


53 


21.28 


-40 


23 


17.7 


35.0 


62.5 


J225322-401931 


22 


53 


22.73 


-40 


19 


31.6 


50.0 


82.0 


J225323-400453 


22 


53 


23.89 


-40 


04 


53.7 


30.0 


10.0 


J225344-401928 


22 


53 


44.89 


-40 


19 


28.6 


167.5 


180.0 


J225404-402226 


22 


54 


04.33 


-40 


22 


26.9 


5.0 


14.0 


J225430-400334 


22 


54 


30.47 


-40 


03 


34.0 


50.0 


8.0 


J225434-401343 


22 


54 


34.70 


-40 


13 


43.2 


5.0 


4.0 


J225436-400531 


22 


54 


36.29 


-40 


05 


31.3 


67.5 


38.0 


J225449-400918 


22 


54 


49.73 


-40 


09 


18.6 


55.0 


66.0 


J225504-400154 


22 


55 


04.67 


-40 


01 


54.8 


5.0 


6.0 



Table 3. Follow-up source list observed in July 2008. 



source name 




R. 


A. 




Dec 


tc/x 


tK 




(12000) 


(12000) 


(min) 


(min) 


J2245 16-40 1807 


22 


45 


16.70 


-40 


18 07.8 


7.0 


4.0 


J224547-400324 


22 


45 


47.88 


-40 


03 24.2 


7.0 


2.0 


J224654-400107 


22 


46 


54.53 


-40 


01 07.2 


7.0 


2.0 


J224750-400148 


22 


47 


50.03 


-40 


01 48.6 


7.0 


2.0 


J224753-400455 


22 


47 


53.73 


-40 


04 55.1 


210.0 


140.0 


J224822-401808 


22 


48 


22.12 


-40 


18 08.2 


7.0 


2.0 


J224827-402515 


22 


48 


27.20 


-40 


25 15.7 


42.0 


4.0 


J2249 19-400037 


22 


49 


19.35 


-40 


00 37.2 


105.0 


28.0 


J224935-400816 


22 


49 


35.21 


-40 


08 16.9 


42.0 


6.0 


J224958-395855 


22 


49 


58.26 


-39 


58 55.4 


14.0 


2.0 


J225004-402412 


22 


50 


04.43 


-40 


24 12.4 


21.0 


6.0 


J225008-400425 


22 


50 


08.82 


-40 


04 25.5 


21.0 


6.0 


J225239-401949 


22 


52 


39.29 


-40 


19 49.5 


42.0 


16.0 


J225529-401101 


22 


55 


29.51 


-40 


11 01.7 


35.0 


8.0 



Table |2] shows the list of 15 sources observed in 2007 with 
target positions and observing times for the C/X band and the 
K-band observations respectively. The 4.8 and 8.6 GHz observa- 
tions were conducted on April 30th, May 1st and 6th 2007. Both 
frequencies were observed simultaneously with a band-width of 
128 MHz. The ATCA operated in the 1.5C aiTay configuration, 
which gives an angular resolution at 5 GHz comparable to the 
one of the original ATESP 1 .4 GHz 6-km ATCA configuration 
(10-15 arcsec). This is important since we want to derive spec- 
tral indices from source flux measurements. Nevertheless we in- 
clude the 6-km antenna in the observations in order to be able 
to perform a higher resolution analysis of the source structure. 
Each source received an average of five linear cuts on a range of 
hour angles between -4^ and +1^ in order to maximize the (m, v) 
plane coverage. A phase calibrator was observed every 20-25 
minutes in 1-min scans. The K-band observations were caiTied 
out on May 25th and 26th 2007 on the same list of sources ob- 
served at 4.8 and 8.6 GHz at two simultaneous IFs (Intermediate 
Frequencies), i.e. at 18.496 and 19.520 GHz. The band-width of 
both IFs was 128 MHz. The total observing time on source is 
reported in Table|2] The array configuration was H214: a hybrid 
and shorter array configuration was chosen in order to map the 
sources with a smaller number of cuts (snapshots) and with an 
angular resolution ( 1 0- 1 5 arc sec) roughly matching the one used 
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Fig. 1. Comparison between ATESP 5 GHz fluxes and new 
5 GHz flux measurements obtained in 2007/2008. Upper limits 
are shown as downwards arrows. Multi-component an/or non- 
Gaussian extended sources are circled. The solid Une represents 
5 (ATESP) = S(2007/2008). 



for the 4.8/8.6 GHz observations. Again the 6-km antenna was 
included in the observations. Each target source was scanned for 
a total of four times during the two-day run in a HA range of 1 ^ 
to4\ 

The list of 14 targets observed in 2008 with the source po- 
sitions and integration times is given in Table [3] The C- and 
X-band run was caiTied out in a single day on July 7th 2008 
with the ATCA in the aiTay configuration 1.5B. An average of 
seven scans was collected for each target source in HA range 
-6^ - -1-4''. The K-band run was scheduled on July 1 1th and 12th 
2008 in two three-hour sessions with the same frequency setting 
and configuration used in 2007. A total number of four cuts in 
the HA range l'^ - 3*^.5 was collected for each target during the 
two observing sessions. The 6-km antenna was included also for 
the 2008 observing runs. 

3.2. Data reduction 

The raw visibility data were reduced using the astronomical 
software package MIRIAD (Sault et al. 1995). The data were 
read into MIRIAD, checked, flagged and split into single-source 
single-frequency files. The bandpass was equalised using the 
standard ATCA calibrator 1921-293 (apart from 2007 K-band 
run where 2243-123 was used instead). The absolute flux scale 
was bootstrapped using the standard ATCA primary calibrator 
1934-638. The phase calibrators 2226-411 and 2232-488 (only 
for the 2008 K-band run) were used to solve for the amplitude 
and phase complex gains in each observing frequency. The cal- 
ibrated visibilities were then imaged by Fourier inversion, af- 
ter excluding the 6-km antenna: the 4.8 and 8.6 GHz visibil- 
ities were imaged separately. The K-band visibilities were in- 
stead combined into a single 19-GHz image by using th e Mult i 
Frequency Synthesis (MFS) technique (Sault & Wieringa ll994l) . 



In the case that the observing frequencies of two or more IFs 
were not too far apart so that radio spectral shape effects were 
still negligible, the MFS technique could be used effectively to 
increase the (m, v) plane coverage and minimise the band-width 
smearing efifects and thus improve the final image quality. As we 
were caiTing out a detection experiment, inverted images were 
inspected using the KARMA package visualisation tool kvis to 
search for at least a 3-sigma detection in the inverted image posi- 
tion were the source was expected to be found (the field centre). 
In case of detection the source images were cleaned and re stored 
using standard clean algorithms: Hog bom (11974 . Clark (Il980l) 
or Steer (Steer, Dewdney, & Ito 1984). The MFS technique was 
used only at the Fourier inversion stage, not at the cleaning stage. 

Limited (m, v) angle coverage of the C/X band observations 
produced elongated restoring beams (24" x 6" at 4.8 GHz and 
14" x3.5" at 8.6 GHz in 2007; 18" x7" at 4.8 GHz and 10" x4" 
at 8.6 GHz in 2008), making the cleaning unreliable. We there- 
fore decided to determine the source fluxes by directly fitting 
a Gaussian point-like model to the calibrated visibilities. In or- 
der to check for the presence of extended structure, we fitted 
the visibilities at different angular resolutions. We compared the 
fluxes obtained by either including or removing the 6-km an- 
tenna, or, for the 8.6 GHz observations only, by removing all 
baselines longer than 750 m (corresponding to 25 k/1) in order to 
have a similar angular resolution as for the 4.8 GHz observations 
(20" x 6"). 

From this analysis we found that most of the sources are 
detected {S IN > 3) at either 4.8 or 8.6 GHz or both, and that 
most of the detected sources can be reliably considered point- 
like at the lowest angular resolution (no 6-km antenna at 4.8 
GHz and no baselines longer than 750 m at 8.6 GHz). For those 
sources (peak) flux densities were determined by fitting a point- 
like Gaussian to the lowest resolution datasets. The remaining 
sources (typically well-resolved or multi-component in the origi- 
nal ATESP 5 GHz catalogue, Prandoni et a/. i2006l) were uv-fitted 
with a single Gaussian whenever suitable, or were alternatively 
fitted with multi-component Gaussians in the restored images 
(convolved to the same resolution as for the other frequencies). 
For such resolved sources integrated flux densities were consid- 
ered to derive the spectral index. 

The (m, v) plane coverage obtained for the K-band observa- 
tions is much better, thanks to the availability of North-South 
baselines of the hybrid ATCA configuration (restoring beam 
~ 10" X 8", excluding the 6-km antenna). Therefore K-band flux 
densities were derived by performing source Gaussian fitting in 
the image plane. Peak (or integrated) flux densities were used for 
unresolved (or resolved) sources to derive the spectral index. 

For all undetected sources 3cr upper limits were derived by 
estimating the noise from V Stokes parameter visibilities for the 
C/X band observations and directly from the images for K-band 
observations. 

In Table |4] we summarise the flux densities at the avail- 
able four frequencies for all the 29 sources observed. The three 
sources not identified with early type galaxies (see discussion in 
Sect. |2]i are explicitely indicated in the Table. For the new ob- 
servations we list also the rms noise values estimated directly 
from the images (for K-band observations) or from V Stokes pa- 
rameter visibilities (for C/X observations). ATESP observations 
were designed to have very uniform image rms values of the or- 
der of 0.08 - 0.09 mJy at 1.4 GH z and 0.06 - 0.08 mJy at 5 
GHz (Prandoni et al. 2000?, Hool). The rms noise values can be 
used to estimate source signal-to-noise ratios and uncertainties in 
the flux density measurements. Table |4] also lists ai and an, the 
spectral indices derived at low and high frequency respectively. 
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Source ATESP Flux 2007/2008 Flux and rms noise values spectral indices 





20 cm 


6 cm 


6 cm 


3 cm 


1.2 cm 








(mJy) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 






J2245 16-40 1807 


9.54 


3.87 


3.30 


0.24 


<1.44 


0.48 


<0.57 


0.19 


-0.71 + 0.04 


<-1.28 


J224547-400324 


32.83 


28.28 


30.49 


0.42 


26.64 


0.60 


27.89 


0.84 


-0.12 + 0.01 


0.06 ± 0.01 


J224654-400107 


5.59 


3.15 


<1.26 


1.00 


<1.44 


0.48 


1.11 


0.26 


-0.45 + 0.05 




J224750-400148 


13.44 


6.09 


5.22 


0.23 


2.63 


0.17 


0.98 


0.30 


-0.62 + 0.02 


-1.25 ±0.29 


J224753-400455 


2.08 


0.67 


1.15 


0.08 


0.43 


0.11 


0.27 


0.06 


-0.89 ±0.19 


-0.59 ± 0.34 


J224822-401808 


19.08 


10.34 


10.04 


0.42 


7.29 


0.61 


4.24 


0.31 


-0.48 + 0.01 


-0.69 ± 0.04 


J224827-402515 


0.58 


0.80 


0.51 


0.17 


1.12 


0.25 


<0.78 


0.26 


0.25 ± 0.27 


<-0.46 


J2249 19-400037 


0.91 


0.64 


<0.33 


0.11 


0.69 


0.16 


0.43 


0.10 


-0.28 ± 0.23 


-0.60 ± 0.32 


J224935-400816 


0.70 


0.82 


1.02 


0.29 


<0.57 


0.19 


<0.60 


0.20 


0.12 ±0.23 


<-0.39 


J224958-395855 


1.52 


1.65 


<0.90 


0.30 


<1.02 


0.34 


1.52 


0.34 


0.06 ±0.12 


_ 


J225004-402412 


3.16 


1.78 


1.89 


0.24 


<0.84 


0.28 


1.28 


0.21 


-0.45 ± 0.08 


-0.28 ± 0.07 


J225008-400425 


2.88 


1.70 


1.61 


0.39 


<0.84 


0.28 


<0.66 


0.22 


-0.41 ± 0.09 


<-0.65 


J225034-401936 


76.62 


25.78 


21.16 


0.30 


17.20 


0.85 


6.03 


0.33 


-0.86 ± 0.01 


-1.33 ±0.02 


J225048-400147('') 


0.60 


0.96 


0.76 


0.17 


1.13 


0.23 


<0.51 


0.17 


0.37 ± 0.27 


<-1.01 


J225057-401522 


2.01 


3.00 


3.44 


0.29 


2.16 


0.40 


2.51 


0.28 


0.31 ±0.08 


0.19 ±0.14 


J225223-401841 


0.98 


0.85 


0.55 


0.12 


<0.51 


0.17 


0.46 


0.13 


-0.11 ±0.20 


-0.13 ±0.21 


J225239-401949 


2.26 


1.18 


1.71 


0.30 


1.02 


0.25 


<0.42 


0.14 


-0.51 ±0.12 


<-1.13 


J225249-401256 


1.52 


1.09 


0.93 


0.15 


0.84 


0.21 


0.71 


0.10 


-0.26 ±0.14 


-0.21 ± 0.24 


J225321-402317r) 


2.32 


1.21 


1.31 


0.11 


<0.57 


0.19 


0.31 


0.09 


-0.51 ±0.11 


-1.05 ±0.15 


J225322-401931 


1.86 


1.01 


0.96 


0.12 


0.92 


0.17 


0.43 


0.09 


-0.48 ±0.14 


-1.06 ±0.23 


J225323-400453 


0.51 


0.85 


0.53 


0.16 


<0.69 


0.23 


0.75 


0.23 


0.40 ±0.31 


0.25 ±0.31 


J225344-401928 


0.60 


3.52 


1.76 


0.10 


2.12 


0.14 


2.05 


0.08 


1.39 ±0.25 


-0.04 ± 0.02 


J225404-402226 


10.34 


3.80 


3.08 


0.43 


3.79 


0.64 


0.78 


0.21 


-0.79 ± 0.03 


-2.01 ± 0.30 


J225430-400334 


<0.26 


0.63 


<0.42 


0.14 


1.43 


0.20 


1.26 


0.29 


>0.70 


-0.16 ±0.21 


J225434-401343 


21.09 


7.80 


4.78 


0.44 


2.95 


0.64 


2.67 


0.43 


-0.78 ± 0.01 


-0.13 ±0.21 


J225436-400531 


0.47 


0.60 


0.90 


0.14 


0.54 


0.15 


<0.45 


0.15 


0.19 ±0.36 


<-0.23 


J225449-400918 


1.24 


0.86 


0.95 


0.15 


0.53 


0.16 


<0.36 


0.12 


-0.29 ±0.17 


<-0.49 


J225504-400154C') 


9.67 


4.59 


4.97 


0.92 


4.91 


1.53 


<1.17 


0.39 


-0.59 ± 0.03 


<-1.82 


J225529-401101 


1.48 


1.09 


1.32 


0.19 


1.34 


0.27 


<0.63 


0.21 


-0.24 ±0.14 


<-0.96 



(") Sources not identified with early-type galaxies. 



In particular we used the original ATESP 1.4 and 5 GHz flux 
densities to derive ai and the new simultaneous 8.6 and 19 GHz 
measurements to derive an- We notice however that for sources 
undetected at 8.6 GHz, was derived by using the simultane- 
ous 4.8 GHz flux measurement, if available. Radio spectra and 
contour plots for the 26 sources associated to early-type galaxies 
are presented at the end of the Paper in Figs. |6] and |2l 

Figure[T]compares the ATESP 5 GHz fluxes (Prandoni et al. 
2006) with the 5 GHz fluxes obtained from the multi-frequency 
observation campaign presented here. The two epoch fluxes 
agree reasonably well with larger scatters for fainter sources 
(generally coiTesponding to lower signal-to-noise ratios). The 
highest deviations from the expected 1 : 1 line at high flux lev- 
els are found for multiple/extended sources (circled dots) where 
new integrated flux determinations, based on poor (m, v)-angle 
coverage data, may be underestimated due to the presence of un- 
detected low surface brightness extended flux. In Sect. |4] we will 
also explore the possible role of source variablity in producing 
significant 5 GHz flux deviations in single-component sources. 

3.3. Full resolution K-band imaging 

ATESP 5 GHz sources were originally imaged at bo th low 
(~ 10") and fun (~ 2") resolution (Prandoni et al. i2006h . 
in order to study their structure (see contours in Fig. |7]i. Here 
we exploit the K-band data to repeat this kind of analysis at a 
higher frequency. K-band images were produced at full resolu- 
tion (~ 2") by including the 6-km antenna in the imaging and 
cleaning process. This has been done for all sources detected in 
K-band low resolution images with 5 /A^ > 6 (7 in total). Results 



Table 5. Full resolution K-band imaging. Source parameters. 



source name 


R.A. 


Dec 




la- 




(J2000) 


(J2000) 


(mJy) 


J224547-400324 


22 45 47.67 


-40 03 23.4 


25.70 


±0.22 


J224822-401808(") 


22 48 21.99 


-40 18 07.8 


3.22 


±0.24 


J225004-402412('') 






1.70 


±0.20 


J225034-401936(*) 






7.15 


±0.31 


J225057-401522 


22 50 57.73 


-40 15 22.3 


2.21 


±0.24 


J225249-401256('') 






0.63 


±0.13 


J225344-401928 


22 53 44.89 


-40 19 27.2 


2.17 


±0.08 


(") Low surface brightness 


flux present to the r 


ight of the source. Including it we get 



S = 4.3 mJy. 

(*) Coordinates are given only for Gaussian fitted sources. Fluxes for multi-component 
and/or low S /N sources were obtained by summing the pixels over the entire source. 



are reported in Table |5] and Fig.|2] where K-band contour plots 
at low (thin contours) and full (thick contours) resolution are 
shown. Unfortunately the full resolution K-band imaging does 
not provide clear morphological information for any source, ex- 
cept for the source J224547 -400324, which shows an elongated 
structure in east-west direction. 



4. Source variability 

Before proceeding with the analysis of the radio spectral prop- 
erties of our sample we should consider the possible effect of a 
source variability in shaping our radio spectra (presented later in 
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22'"45"49?0 48=5 4B?0 47?5 47?0 asHs^EafO 33=5 33fO 31=5 2lf0 33''50™5?5 5f0 4=5 4=0 3?5 SS^SO^SS^S 35'0 34=5 34=0 33=5 

RA (J300D) RA (J3000) RA (J3000) RA (J3000) 




33''50"'59=0 53=5 58=0 57=5 57=0 50=5 33''53°'51=0 50=5 50=0 49=5 49=0 33^53"46=0 45=5 45=0 44=5 44=0 

R.A (J2000] RA (jaOOO) RA [J2000) 



Fig. 2. K-band full resolution contours (thick lines) superimposed on low resolution contours (thin lines) for the seven sources 
detected in low resolution K-band images with S /N > 6. Contour levels are 3, 4, 5, 6, 8, 10, 15, 20, 30, 50 and 60 % of the peak 
flux density. 



Fig.|6l). Although high frequency observations were done simul- 
taneously in either 2007 or 2008, low frequency observations 
were done long before that, between 1999 and 2005. 

From Sadler et al. (|2006|) we know that high frequency (K- 
band) variability is not very important on time scales of a few 
years: only 5 % of the sources turn out to be variable at a level 
> 30 % over a baseline of 1-2 years. We assume that at the lower 
frequencies (1.4 and 4.8 GHz) the variability will be lower than 
that, since the low frequency emission may mostly come from 
extended, optically thin regions. However, we cannot exclude 
that in ten years time the overall level of emission at high fre- 
quencies has changed enough so that, for example, the spectral 
index qm (between 5 and 8.6 GHz) may have varied while 
(between 1.4 and 5 GHz) and oh (between 8.6 and 19 GHz) re- 
mained (roughly) constant. Such a change may manifest itself 
if, for example, both the high and low frequency spectral indices 
are negative, but the intermediate spectral index is positive. This 
kind of effect does indeed happen, as a visual inspection of Fig.|6] 
easily reveals (see e.g. sources J225057-401522 and J225404- 
402226), although it is not so big as to influence the discussion 
in a significant way. This co nfirms and strengthens the conclu- 
sion of Sadler et al. (2006) that strong source variability is a 
minor effect, even on time scales of the order of 5-10 years. 

For a more quantitative analysis of the source variability in 
our sample, we compared the original ATESP 5 GHz fluxes (ob- 
tained in 2000/2001, see Prandoni et al. 2006) with the ones 
obtained with the 2007/2008 follow-up and checked the signif- 
icance of their variations through the analysis of the so-called 



variability index. In particular we analysed two well-know n vari- 
ability indices, defined as follows (see Gregorini et al. 119861 : 
Quirrenbach et al. 1 19921) : 

Vi = AS/cTto, 
V2 = AS/ <S > 

where AS = S i - S2, with Si and 5 2 being the 5 GHz fluxes 
measured respectively in epochs 1 (2000/01) and 2 (2007/08); 
a-,„, = ^JcT^{Sl) + o-\S2), and < 5 We notice that 

V2 represents the so-called modulation index as it reduces in the 
case of two epochs measurements. 

Figure |3] shows the Vi and V2 variability index distribution 
as a function of the signal-to-noise ratio for single-component 
sources (multi-component sources are unlikely to vary). To com- 
pute Vi we estimated cr(5 1) and cr(52) through Condon's master 
equations (Condon 1997) by assuming calibration errors of the 
order of 3% for both epochs (see Prandoni et al. .2000^) . Dotted 
and solid lines in Fig.[3]correspond to random Gaussian AS fluc- 
tuation probabilities of 0.5% and 0.1% respectively. In general 
a source is considered as a good variable candidate whenever 
the above probability is < 0.1% (or AS > 3>.29cr,oi)- It is clear 
from the figure that two sources (J224753-400455 and J225344- 
401928) satisfy both the V\ and V2 variability criteria, while a 
third source (J224547 -400324) is defined as variable only with 
respect to V2 and for this reason has to be considered a less reli- 
able candidate. All the other sources have 5 GHz flux variations 
consistent with being random fluctuations. 
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Fig. 3. Variability index Vi (Top) and V2 (Bottom) for single- 
component sources as a function of S/N (referred to 2007/8 ob- 
servations). Arrows indicate sources which were not detected in 
epoch 2. Dotted and solid lines correspond to random AS fluc- 
tuation probabilities of 0.5% and 0.1% respectively. 



From the analysis of the overall source spectra (see Fig. |6]l 
we see that source J224753-400455 (AS = -4.7cr,o,; AS = 
-0.53 < S >) has a smooth power-law spectrum from 1.4 
to 19 GHz, with the only exception of the 5 GHz point ob- 
tained in 2008. This argues against variability for this sources. 
On the other hand source J225344-401928 (AS = ll.lcr,^,; 
AS = 0.67 < S >), clearly shows a change in the spectral 
shape from strongly inverted to flat when we go from the ATESP 
non-simultaneous flux measurements (red squares) to the present 
quasi-simultaneous flux determinations (blue diamonds). This 
strongly supports the hypothesis of variability. Source J224547- 
400324 is a less reliable candidate (AS = -1 .7cr,£„) and flux vari- 
ations, if real, are quite small (AS = -0.08 < S >). The overall 
spectrum is rather flat over the entire frequency range covered. 
Both the flatness of the spectrum and the radio morphology (not 
inconsistent with a core-jet structure, see Fig.|2|i may support the 
presence of some level of variability in this source. 

In summary, we find that one of the 21 (e.g. 5%) single- 
component sources in our sample is a strong variable, while 
a second source may present some low-level var iabilit y. This 
is consistent with previous results. Tingay et al (l2003b found 
that the level of variability of powerful (~ 1 Jy), compact ra- 
dio sources moderately increases with frequency from a median 
variability of 6% at 1.4 GHz to 9 % at 8 .6 GHz over a time-scale 
of 3 - 4 years. Barvainis et al. (|2005) found similar variabil- 
ity levels at 8.5 GHz for a sample of radio-loud and radio-quiet 
QSOs observed with the Very Large Array. It is worth noting 
that the flux density levels probed by the radio-quiet QSOs of 
the Barvainis et al. sample are similar to the flux densities of our 
sources (; 1 mJy). 



5. Discussion of radio spectra 

Figure |6] shows the overall radio spectra of the 26 ATESP 5 GHz 
sources associated with early-type galaxies discussed in this 
work. Red circles represent non-simultaneous ATESP 1 .4 and 5 
GHz observations, while blue diamonds represent the new quasi- 
simultaneous observations at 4.8, 8.6 and 19 GHz. We notice 
that low frequency spectral behaviour (steep or flat) derived from 
1 .4 and 5 GHz ATESP data is generally confirmed by high fre- 
quency simultaneous observations. This can be considered as a 
first qualitative confirmation of the fact that the sub-mJy flat- 
/inverted-spectrum population associated with early-type galax- 
ies revealed by the ATE SP 1.4 and 5 GHz surveys (Prandoni et 
al. I2OO6I: Mignano et al. '2008) and discussed in Sect. [T] is gen- 
uine. The most notable exception is represented by the variable 
source J225344-401928 (see Sect.|4|, which is remarkably char- 
acterised by a very inverted spectrum between 1.4 and 5 GHz 
(ai = 1 .4), not confirmed as real by high frequency simultane- 
ous observations (om ~ uh ~ 0). 

For a more detailed analysis of the radio spectra we will fo- 
cus below on two spectral indices, ai and uh (listed in Table|4]i. 
These two spectral indices can be used in a sort of "two-colour" 
diagram as was done in previuos works (see Sadler et al. 200^ 
Massardi et al 2008; Tucci et al. 2008). Note that the analysis of 
the two-colour diagram is not exhaustive, and the spectral index 
between 4.8 and 8.6 GHz (e.g. um) may give further independent 
information on the spectrum, as previously noted in Sect. |4] 

We will also compare our data with recent high frequency 
sur veys (Sadler et al. 12006; Massardi et al. 12008 ; Tucci et 
al. I2OO8I) and will try to determine how and to what extent 
frequency selection and flux limits influence the results. The 
ATESP 5 GHz survey is the only one selected at low frequency, 
while the other three are based on 20 - 30 GHz observations. 
Also there is quite a difference in the flux limit, which reaches 
from ~ 0.6 mJy (our sample) via 20 mJ y at 33 GHz (Tucci et al. 
2008), 100 mJy at 20 GHz (Sadler et al. |2006|) to 500 mJy at 20 
GHz (Massardi et al. I2008D . One extra selection criterion used 
for our sample was that we followed-up at high frequencies only 
sources associated with galaxies having early-type spectra; for 
that reason we also defined a sub-sample of 36 sources associ- 
ated with galaxies from the Bright Source Sample discussed by 
Massardi et al. (2008). Since the latter sources are all bright in 
the radio, these are presumably all classical radio galaxies; their 
typical radio powers put most of them effectively in the FRII 
range. 

Figure |4] shows two-colour diagrams of our (red symbols) 
and the Massardi et al. galaxy (blue diamonds) samples. Also 
plotted are all sources from the Massardi et al. sample that are 
associatd with quasars or quasar candidates (blue dots). We di- 
vided the sources in a nearby (z < 0.5, left panel) and a far 
(z > 0.5, right panel) sample, in order to check if distance may 
have an effect upon the spectral indices. In the far sample we 
included quasar candidates with unknown redshift, the large ma- 
jority of which (presumably) are at redshifts > 0.5 on the basis 
of their optical magnitude distribution. 

For a physical interpretation of the two-colour diagram we 
divided it in four quadrants, as done in previous works: 1) 
sources with ai < 0, an < 0, conventionally named steejQ, 
2) sources with ai < 0, an > 0, named upturn; 3) sources 
with ai > 0, an > 0, named inverted; 4) sources with ai > 
0, an < 0, named peaked. We point out that the radio spectral 
index of the overall synchrotron emission from a jet-dominated 

' In fact in this quadrant fall both steep (a < -0.5) and flat (-0.5 < 
a <0) sources 



8 



I. Prandoni et al: The ATESP 5 GHz radio survey. Ill 

z<0.5 z>0.5 / no z 





Fig. 4. High- against low-frequency spectral index. Red symbols indicate ATESP sources: filled symbols are used for an derived 
between 8.6 and 19 GHz and empty symbols are used for an derived between 4.8 and 19 GHz; triangles indicate upper/lower limits 
in an/ai respectively. Blue diamonds indicate sources associated with early type optical spectra from the M assardi et al. (20081) 
sample; blue dots indicate sources associated with quasars or quasar candidates in the Massardi et al. (l2008l) sample. The typical 
error box for the Massardi et al. sources is indicated in the upper left corner. Left panel: objects with redshift < 0.5. Right panel: 
objects with either redshift > 0.5, or unknown redshift. The dashed line represents the location of sources with power-law spectra 
(i.e. with Oh - ai)- Horizontal and vertical dotted hnes respectively indicate the an -0 and = lines defining the four quadrants 
discussed in the text. 



radio source is expected to be from steep to moderately inverted 
(-0.7 < a < 0.2) at the frequencies of this work, depending on 
the relative contributions of extended (optically-thin) and base 
(self-absorbed) jet components. Jet-dominated sources are there- 
fore expected to mostly fall into the steep quadrant. Instead, for 
an ADAF, one exp ects 0.2 < a < 1.1 up to mm wavelengths (see 
Nagar et a/. 1200 lb . with a varying according to the accretion rate 
(f.i. a = 0.4 when L ~ 10""*LEdd and a ~ 1 when L < lO^^Lfidd; 
Mahadevan ll997h . ADAF systems are therefore expected to fall 
into the inverted quadrant. Co-existence of moderate outflows 
can flatten the spectrum (with a remaining positive), whereas 
strong outflows can shift the peak of the radio emission from 
mm to cm wavelengths (Quataert & Narayan 1999). 

As can be seen in Fig. |4] our ATESP sources, as well as 
the sources from the much brighter Massardi et al. galaxy sam- 
ple, are mostly located in the lower part of the plots: 23/24 and 
33/36 sources for the ATESP and the Massardi sample respec- 
tively have Oh < 0.2. Although these numbers are significant 
only at about the 3cr level due to the relatively poor statistics in- 
volved, there can be little doubt that both samples consistently 
contain more sources with an < 0.2. This is of course hardly 
surprising for the ATESP sample, which was originally selected 
at 5 GHz, but the high frequency selection of the Massardi sam- 
ple does not make much difference, at least in this respect. Less 
obvious a priori is the general steepening (ai > an) of the spec- 
tra from low to high frequency. This is evident from the large 
number of sources lying below the diagonal dashed line, cor- 
responding to power-law spectra (19/24 and 31/36 for ATESP 
and Massardi). These numbers are significant only at the 2cr 
level, but again both samples behave consistently and actually 
appear to cover identical regions in the two-colour diagram. 
Both facts favour a jet-dominated scenario, with steep-spectrum 
{an < -0.5) sources dominated by synchrotron optically-thin 



Table 6. Mean values and standard deviations of the low and 
high frequency spectral indices. 



Sample*"' cr cr^ a 

ATESP ^OTS (OTT) =079 ((TTeT 

Massardi G -0.01 (0.08) -0.35 (0.09) 

ATESP + Massardi G far -0.08 (0.07) -0.50 (0.08) 

Massardi Q far 0.21 (0.04) -0.16 (0.03) 

'Massardi G' and 'Massardi Q' respectively refer to galaxies and 
stellar identifications of the Massardi et al. sample. 'Far' refers to 
z > 0.5 sources (or unknown redshift for quasar candidates) 



emission and flat-spectrum {an > -0.5) sources dominated by 
synchrotron optically-thick emission coming from the base of 
the jet. Nevertheless jet+ADAF models can also be consistent 
with values Qlan'O.!. Pure ADAF models {ai and an > 0.2) 
are ruled out by the present data (only one object in the Massardi 
galaxy sample can be considered as a pure ADAF candidate), in- 
dicating that such radiatively ineflicient accretion regimes must 
be very rare also at the sub-mJy levels probed by the ATESP 

5 GHz sample. 

For a more quantitative analysis of the two-colour diagram 
we statistically compared both the low- and high-frequency 
spectral index distribution of the ATESP and the Massardi et al. 
galaxy samples by implementing the survival analysis methods 
described in Feigelson & Nelson (il985h and Isobe et al. (Il986h . 
This allows us to properly treat upper/lower limits in anlai re- 
spectively. Both ol and qh spectral index distributions are con- 
sistent with being drawn from the same parent population. In fact 
this hypothesis is rejected at only a < 2cr confidence level (Peto 

6 Prentice generalized Wilcoxon test, Isobe & Feigelson 1990i). 
This is especially interesting when we consider that the ATESP 
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and the Massardi et al. samples cover very different ranges in 
flux density (by two or three orders of magnitude). 

A similar comparison can be made between galaxies and 
quasars. To improve our statistics we consider the ATESP and 
the Massardi et al. galaxies as a single population, and we 
limit our analysis to z > 0.5 sources, where most quasars (or 
quasar candidates) are located (see right panel of Fig. HJi. In 
this case we find that the hypothesis that z > 0.5 galaxies and 
quasars are extracted from the same parent population is re- 
jected at a > 4.5cr level (Peto & Prentice generalized Wilcoxon 
test, Isobe & Feigelson 119901) . Table summarises the mean 
values for the samples discussed above (Kaplan-Meier estima- 
tor, Isobe & Feigelson 1990). We notice that the Massardi sub- 
sample of z > 0.5 quasars and quasar candidates (Massardi Q 
far) tend to be flatter or even more inverted than sources associ- 
ated with galaxies, indicating that quasar-like objects are more 
commonly associated with compact sources and less with the 
classical extended radio lobes in which the optically thin syn- 
chrotron radiation with steep radio spectra is dominant. On the 
other hand, no difference is seen in spectral steepening: if we 
compare the curvature C = an - ai (see e.g. Gregorini et al. 
11984 . we find C = -0.42 + 0.11 (ATESP + Massai-di G far), 
and C = -0.37 + 0.05 (Massai'di Q far). 

Following Sadler et al. ('2006') we consider in more detail the 
distribution of sources in the four quadrants of the two-colour di- 
agram, which represent different types of radio sources (defined 
earlier in this section). The distribution of objects over the four 
quadrants is shown in Table |7] (where U stands for Upturn, I for 
Inverted, S for Steep , P for Peaked). Most ATESP sources are in 
the lower left quadrant (61%, Table [T), as for classical FRI and 
FRII radio galaxies. This may be due to the early type galaxy 
optical pre-selection, and in fact the entries of the Massardi et 
al. galaxy sample (M-G in Table |7ll are, considering the errors, 
practically the same as those of the ATESP sample, despite the 
very different flux range covered by the two samples. 

In Table [T] we also compare the distribution of objects in 
three samples, all selected at high (> 20 GHz) frequencies, 
where no optical pre-selection was applied (Tucci et al. l2008t 
Sadler et al. 2006; Massardi et al. 2008). There appear to be 
some trends with flux density, which we can attribute partly 
to selection effects and partly to the different composition of 
the overall samples, in which quasar spectra may be more or 
less dominant, depending on the limiting flux density. The up- 
turn and steep spectra (i.e. the left hand side of the two-colour 
diagram) tend to become more numerous at lower flux densi- 
ties, while the inverse happens with inverted and peaked spectra 
(right hand side). These trends are individually significant only 
at the 2cr level, again due to the small numbers involved, but ap- 
pear to be quite consistent in all four quadrants. If we sum the 
U and S sources (those with ai < 0) and the / and P sources 
(ai > 0), the dependence on the flux density becomes quite ob- 
vious: whereas about 18 % of the sources have a positive (I 
and P) at the flux density level of the Tucci et al. (2008) sam- 
ple, at the much higher flux levels of the Massardi et al. sample 
this strongly increases to about 60 %; this difference is definitely 
significant at about Act. The fact that the number of rising spec- 
tra inc reases with flux limit was already noticed by Bolton et al. 
( |2004|) . No doubt these trends may be due to selection effects; 
however, the possible increase of peaked sources at higher flux 
densities is harder to explain by selection effects only. 



Table 7. Spectral type statistics in five samples' 





ATESP M-G 


Tuc 


Sad 


M 




> 0.6 mJy > 500 mJy 


> 


20 mJy 


> 100 mJy 


> 500 mJy 






(percent) 








(percent) 




u 


9 


±5 8 


±5 


29 


±6 


22 ±5 


10 +1 


I 


13 


7 12 


6 


7 


3 


20 4 


28 5 


s 


61 


14 47 


11 


53 


7 


34 6 


30 5 


p 


17 


8 33 


11 


11 


3 


24 5 


32 5 



ATESP 5 GHz sample (26 sources); Massardi et al. galaxy 
sample (M-G, 36 sources); Tucci et al. (2008) sample (Tuc, 
102 sources); Sadler efaZ. (2006) sample (Sad, 101 sources) 
Massardi et al. (2008) overall sample (M, 218 sources) 



6. The AGN component of the ATESP 5 GHz sample 

In Sect. |5]we discussed the radio spectral properties of a repre- 
sentative sub-sample of ATESP radio sources associated to early 
type galaxies in comparison with some recent high-frequency 
radio surveys. The main results of such an analysis can be sum- 
marised as follows. Most sources have spectra ranging from 
steep to moderately inverted, which is consistent with them be- 
ing jet-dominated systems. Strong spectral similarities found be- 
tween our and the much brighter Massardi et al. galaxy sam- 
ple may suggest that our radio sources are simply lower lu- 
minosity counterparts of FRII radio galaxies. No evidence for 
an ADAF population was found, and peaked GPS-like sources 
seem to be preferentially found at higher flux densities than 
the ones probed by the ATESP. In addition quasar-like upturn 
sources are rare (9 + 5%, see Table |7J. While this latter result 
is partly a consequence of the ATESP optical pre-selection, it 
argues against an obscured radio-quiet quasar population hid- 
den among the ATESP early type galaxies. Nevertheless com- 
pact quasar-like objects may be present among flat/moderately 
inverted jet-dominated sources. 

To further explore this hypothesis and probe the AGN com- 
ponent of the ATESP sample at large we now compare linear 
sizes and radio powers of all the ATESP sources with redshift de- 
termination (71 source s, see Table [T]i, wi th the ones of the well- 
know n B2 (CoUa et al. 119751: Fanti et al. 1 1978b and 3CR (Laing 
et al. 11983 ') radio source catalogues, and with sizes and powers 
of a s ample of optically-selected radio-quiet QSO (Leipski et al. 
I2006h . B2 and 3C catalo gues probe different ranges in radio lu- 
minosity, with the former dominated by FRI (P < 10^^ W/Hz) 
and the latter by FRII (P > 10^^ W/Hz) radio galaxies. Such 
a comparison is shown in Fig. |5] (left panel). In the right-hand 
panel of Fig.|5]the same plot is shown, where symbols now stand 
for different optical types of ATESP radio sources. 

We first focus our analysis on ATESP resolved sources (red 
squares in left panel). Such sources have generally larger sizes 
than radio-quiet QSOs and seem to be more consistent with the 
B2 radio galaxy size distribution. In particular they seem to con- 
firm and extend to lower radio luminosity the size-power relation 
found for B2 radio galaxies (de Ruiter et aZ. ,1990) . When enlarg- 
ing the analysis to unresolved sources (red upper limits in left 
panel), we notice a group of sources at f > 10^^ W/Hz which 
clearly stand below the B2/3C source size distribution. As shown 
in the right panel of Fig. |5] such sources are in fact associated 
to broad-line AGNs (see green upper limits) and could be gen- 
uine compact sources, perhaps the radio-intermediate counter- 
parts of radio-quiet QSOs. Unresolved sources with radio pow- 
ers P < 10^^ W/Hz are mostly associated to early-type galaxies. 
In general the given upper limits do not allow us to distinguish 
between B2 FRI-like and quasar-like sizes. Nevertheless in some 
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Fig. 5. Linear size as a function of radio power (at 1 .4 GHz). Left panel: comparison between four different samples: ATESP 5 GHz 
sample (red squares and upper limits), B2 catalogue (blue open circles), 3C catalogue (black filled circles), radio-quiet QSO (green 
asterisks and upper limits). For Leipski et al. radio-quiet QSO sample 1.4 GHz luminosities were derived from tabulated 5 GHz 
values, assuming a = -0.7. Right panel: comparison between different types of objects in the ATESP 5 GHz sample: sources 
associated to early type galaxies (red squares and upper limits); sources associated to late type galaxies (blue open circles and upper 
limits; among late type sources we may also have some mis-classified naiTow-line AGN); sources associated to broad-line AGNs, 
typically quasars (green asterisks and upper limits). 



cases we have dllO kpc, making such sources more consistent 
with a radio-quiet scenario. Some of those sources may also be 
low power BL Lac, where the optical activity is very weak (see 
discussion in Sect. [Hi. From geometrical considerations we ex- 
pect < 3-6% (or < 2-4) of such objects in the samp le, assuming 
viewing angles < 15°-20°(Zenzus & Pearson [l987|). Such frac- 
tion could be somewhat higher, due to possible sample contam- 
ination from boosted objects belonging to a parent population 
fainter than the ATESP flux limit. For a more conclusive analy- 
sis of the size-power plot we need higher resolution data for the 
ATESP sample. In particular VLBl imaging would allow us to 
better probe the low power BL Lac scenario. As a final remark it 
is worth noticing that ATESP source sizes may suffer from some 
underestimation. This because low surface brightness flux asso- 
ciated to radio jets could be missed by blind source Gaussian 
fitting procedures performe d to construct the ATESP source cat- 
alogues (see Prandoni et Q/. l2000btl2006l) . 

7. Summary 

We presented follow up quasi-simultaneous observations at 4.8, 
8.6 and 19 GHz for a representative sub-sample of ATESP 
sources associated to early type galaxies (26 sources with S > 
0.6 mJy). The new data were discussed together with the original 
1.4 and 5 GHz ATESP data (Prandoni et aZ. l2000bi: l2006l) . with 
the aim of better understanding the radio spectral properties of 
such sources and ultimately the nature of AGNs at sub-mJy flux 
levels. In particular we were interested in understanding whether 
the AGN component of the sub-mJy population is more related 
to efficiently accreting systems - like radio-intermediate/quiet 
quasars - or to systems with low accretion rates - like e.g. FRI 



radio galaxies - or to low radiative efficiency accretion flows - 
like e.g. ADAF. 

From the analysis of the radio spectra, we found that our 
sources are most probably jet-dominated systems, with steeper 
sources dominated by optically-thin synchrotron emission (typ- 
ical of FRI and FRII radio galaxies) and flatter or moderately- 
inverted sources dominated by optically-thick emission coming 
from the base of the jet. Pure ADAF models seem to be ruled out 
by the high frequency data, while ADAF+jet scenarios could be 
consistent with flat/moderately inverted sources, but are not re- 
quired to explain the data. 

By comparing our sample with some recent high frequency 
(;20 GHz) surveys we found strong spectral similarities be- 
tween our and the much brighter Massardi et al. galaxy sample, 
possibly suggesting that our radio sources are mostly lower lu- 
minosity counterparts of the Massardi et al. (FRII) radio galax- 
ies. In addition peaked GPS -like sources seem to be prefer- 
entially found at higher flux densities than the ones probed 
by the ATESP sample, and quasar-like upturn sources are rare 
(9 + 5%). This latter result is partly due to the ATESP early type 
galaxy pre-selection, but argues against an obscured radio-quiet 
quasar population hidden among the ATESP early type galaxies. 
Nevertheless compact quasar-like objects may be present among 
flat/moderately inverted jet-dominated sources. 

To further investigate such a hypothesis we compared radio 
powers and linear sizes of the ATESP 5 GHz sample with the 
ones of two well-known brighter radio source catalogues (B2 
and 3CR) and with the sample of radio-quiet quasars studied 
by Leipski et al. (2006). Unfortunately such a comparison suf- 
fers from the large number of unresolved ATESP sources (or 
the large number of size upper limits). A significant fraction of 
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Fig. 6. Flux densities (in mJy) against frequency (in GHz) for the 26 ATESP 5 GHz sources associated to early-ty pe gala xies 
discussed in t his wo rk (see Table [1]). Red filled circles represent the original 1.4 and 5 GHz points (Prandoni et al. 2000b and 
Prandoni et al. |2006|) . while blue diamonds represent the 2007 and 2008 observations discussed in the present paper. No error bars 
are drawn, since these are typically smaller than symbols. Dotted lines indicate the expected slopes for a = and a - -0.5. 



ATESP sources have sizes consistent with the ones of B2 and 
3CR radio galaxies. Such sources seem to confirm and extend 
to lower radio luminosity the size-p ower relation found for B2 
radio galaxies (de Ruiter et al. Il990l) . On the other hand, a com- 
ponent of very compact d < 10 kpc radio sources associated to 
early type galaxies seem to have sizes more consistent with the 
ones of radio-quiet quasars. Among them we may have a few low 
power BL Lac. For a more conclusive analysis we need higher 
resolution (possibly VLBI) data for the ATESP sample. 



Acknowledgements. The Australia Telescope Compact AiTay is part of the 
Australia Telescope which is funded by the Commonwealth of Australia for op- 
eration as a National Facility managed by CSIRO. 



References 

Abramowicz, M.A.. Igumenshchev, I.V., Quataert, E., & Nar ayan, R. 2002, 
ApJ, 565, 1101 

Afonso, J., Mobasher, B., Koekemoer, A., Norris, R. R, & Cram, L. 2006, AJ, 
131, 1216 

Barvainis, R., Lehar, J., Birkinshaw, M., Falcke, H., & Blundell, K. 2005, ApJ, 
618, 122 



12 



I. Prandoni et al: The ATESP 5 GHz radio survey. Ill 



1-3 



X 



a. 



100 
10 

1 

0.1 
100 

10 

1 

0.1 
100 

10 

1 

0.1 
100 

10 

1 

0.1 



, , t , 

3^ — I I I I llll| — H 



-| 1 — I — I I I 1 1 1 1 1 — r 

J225239-401949 



■-0----;. 



J225323-400453 



J225430-400334 



-1—1 — 1 1 1 1 



o o 



!1 M I I I M il 



J225449-400918 



J I I 'A L_ 



-| 1 — I I I I 1 1 1 1 1 — r 

J225249-401256 



+m 



^ — I 1 1 1 iiii| — h 



J225344-401928 



^\ — I I 1 1 lllll — \—\-^\ — I I 1 1 lllll — \—\-^ — I I 1 1 lllll — H 



0""0 



J225434-401343 



I I I I llll| \—^ 

J225529-401101 



J I I I L 



T — I — I I I 1 1 1 1 1 — r 

J225322-401931 



-<»-0:. 



■o 



J225404-402226 



O 



J225436-400531 



wL 



I I I I 1 1 1 1 1 I I 



frequency (GHz) 



Fig. 6. Continued. Flux densities (in mJy) against frequency (in GHz) for the 26 ATESP 5 GHz sources associated to early-type 
galaxies discussed in this work (see TablelTJ- 



Bolton, R.C., Cotter, G., Pooley, G.G., et al. 2004, MNRAS , 354, 485 

Clark, B.G. 1980, A&A, 89, 377 

Colla, C, Fand. C, Fanti, R., et al. 1975, A&AS, 20, 1 

Condon, J.J. 1984, ApJ, 287, 461 

Condon, J.J. 1997, PASP, 109, 166 

De Zotti, G., Ricci, R., Mesa, D., et al. 2005, A&A, 431, 893 

Doi, A., Kameno, S., Kohno, K., Nakanishi, K., & Inoue, M. 2005, MNRAS, 

363, 692 

Falcke, H. & Biermann, PL. 1999, A&A, 342, 49 

Fanaroff, B. L. & Riley. J. M. 1974, MNRAS, 167, 31 

Fanti, R., Gioia, 1., Lari, C, Ulrich, M.H. 1978, A&AS, 34, 341 

Feigelson, E.D. & Nelson, PI. 1985, ApJ, 293, 192 

Georgakakis, A., Mobasher, B., Cram, L., et al. 1999, MNRAS, 306, 708 

Giroletti, M., Giovannini, G., & Taylor, G.B. 2005, A&A, 441, 89 

Gregorini, L., Mantovani, R, Eckart, A., et al. 1984, AJ, 89, 323 

Gregorini, L., Ficarra, A., & Padrielli, L. 1986, A&A, 168, 25 

Gruppioni, C, Ciliegi, P, Rowan-Robinson, M., et al. 199 9, MNRAS, 305, 

297 

Isobe, T. & Feigelson, E.D. 1990, BAAS, 22, 917 
Isobe, T, Feigelson, E.D. & Nelson, PI. 1986, ApJ, 306, 490 
Jarvis, M.J. & Rawlings, S. 2004, New Astron. Rev., 48, 1 173 
Hildebrandt, H., Erben, T, Dietrich, J.P, et al. 2006, A &A, 452, 1121 
Hogbom, J. A. 1974, A&AS, 15, 417 



Klockner, H.-R., Martfnez-Sansigre, A., Rawlings, S., & Gan'ett, M.A., 2009, 
MNRAS, in press 

Lacy, M., Laurent-Muehleisen, S.A., Ridgway, S.E., Becker, R.H.. & White, 
R.L. 2001, ApJ, 551, L17 

Laing, R.A., Riley, J.M., & Longair, M.S. 1983, MNRAS, 204, 151 

Leipski, C, Falcke, H., Bennert, N., & Huttemeister, S. 2006, A&A, 455, 161 

Magliocchetti, M.. Maddox, S. J., Wall, J. V., Benn, C. R. , & Cotter. G. 2000, 

MNRAS, 318, 1047 

Mahadevan, R. 1997, ApJ, 477, 585 

Massardi, M., Ekers, R.D., Murphy, T, et al. 2008, MNRAS, 384 , 775 
Mignano, A., Prandoni, I., Gregorini, L., et al. 2008, A& A, 477, 459 
Muxlow, T.W.B., Richards, A.M.S., Garrington S.T., etal. 2005, MNRAS, 358, 
1159 

Nagar, N.M., Wilson, A.S., & Falcke, H. 2001, ApJ. 559. L87 
Narayan, R. & Yi, I. 1994, ApJ, 428, L13 

Nonino, M., Berlin, E.. da Costa, L., et al. 1999, A&AS, 137, 51 
O'Dea, C.R 1998, PASP, 1 10, 493 

Olsen, L.F., da Costa, L., Madejsky, R., et al. 2006, A&A , 456, 881 
Prandoni, I., Gregorini, L., Parma, P, etal. 2000, A& AS, 146, 31 
Prandoni, I., Gregorini, L., Parma, P, etal. 2000, A& AS, 146, 41 
Prandoni, I., Gregorini, L., Parma, P, etal. 2001, A& A, 365, 392 
Prandoni, I., Gregorini, L., Parma, P, etal. 2001, A& A, 369, 787 
Prandoni, I., Parma, P, Wieringa, M.H., et al. 2006, A&A , 457, 517 
Quataert, E. & Narayan, R. 1999, ApJ, 520, 298 



I. Prandoni et al. : The ATESP 5 GHz radio survey. Ill 



13 



Quirrenbach, A., Witzel, A., Kirchbaum, T.P., et al. 1992, A&A, 258, 279 

Ricci, R., Sadler, E.M., Ekers, R.D., et al. 2004, MNRAS, 354, 305 

Richards E.A., Fomalont, E.B., Kellermann, K.I., et al. 1 999, ApJ, 526, L73 

de Ruiter, H.R., Parma, P., Eanti, C, & Fanti, R. 1990, A&A, 227, 351 

Sadler, E.M., Ricci, R., Ekers, R.D., et al. 2006, MNRAS, 371, 898 

Sault, RJ. & Wieringa, M.H. 1994, A&AS, 108, 585 

Smolcic, v., Schinnerer, E., Scodeggio, M., et al. 2008, ApJS, 177, 14 

Snellen, I.A.G., Schilizzi, R.T., Miley, G.K., et al. 2000, MNRAS, 319, 445 

Steer, D.G., Dewdney, PE., & Ito, M.R. 1984, A&A, 137, 159 

Tingay, S.J., Jauncey, D.L., King, E.A., et al. 2003, PASJ, 55, 351 

Tucci, M., Rubino-Martin, J.A., Rebolo, R., et al. 2008, MN RAS, 386, 1729 

Vettolani, G., Zucca, E., Merighi, R., et al. 1998, A&AS, 130, 323 

Waldram, E.M., Bolton, R.C., Pooley, G.G., & Riley, J.M. 2007, MNRAS, 379, 

1442 

Zenzus, J.A. & Pearson T.J. 1987, in 'Superluminal Radio Sources', 
Cambridge University Press, Cambridge, p. 1 



14 



I. Prandoni et al. : The ATESP 5 GHz radio survey. Ill 




Fig. 7. Contour images of the 26 ATESP 5 GHz sources identified with early-type galaxies, discussed in this work. Thin black con- 
tours represent low resolution 5 GHz data; superimposed are high resolution 5 GHz contours (thick black lines) and low resolution 
K-band contours (thick gray fines). Contour levels are 2.5 (only K-band), 3, 4.5, 6, 10, 20, 50 and 100 % of the peak flux density. 
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Fig. 7. Continued. Contour images of the 26 ATESP 5 GHz sources identified with early-type galaxies, discussed in this work. 



